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I. INTRODUCTION 

We examine applicability of Schwinger-Dyson equa- 
tions in simple models characterized by a nontrivial 
vacuum. An infinite set of Schwinger-Dyson equations 
(SDE's) represents integral relations between Green's 
functions that in principle describe the complete dynam- 
ics of the underlying field theory. In terms of loop ex- 
pansion even a single SDE contains an infinite series of 
interaction terms. For this reason, in QCD, where strong 
interactions between quarks and gluons dominate long 
range dynamics Schwinger-Dyson equations have been 
extensively used to describe various non-perturbative 
phenomena; ranging from confinement and chiral sym- 
metry breaking, to applications in hadron phenomenol- 
ogy [Il-O]- Even when the underlying theory has only 
a limited number of elementary interactions the full set 
of SDE's generates a complicated effective potential. In 
practical applications any approximation to SDE's elim- 
inates an infinite set of such effective interactions and 
therefore it is important to access the applicability of 
any such truncation in QCD phenomenology. A num- 
ber of investigations in the ultraviolet, and in the more 
relevant for strong QCD, the infrared region, have been 
performed [Hll]. 

There is ample evidence from lattice simulations that 
confinement in the QCD vacuum has origin in topol- 
ogy 6-9 . It has been postulated long ago that both con- 
finement and chiral symmetry breaking originate from in- 
stantons. In the case of confinement topological objects 
like center vortices or magnetic monopole loops perco- 
late through Wilson loops and lead to its area law de- 
pendence [Tntil3j . As such a condensate of magnetic 
monopoles ought to screen electric flux lines and pro- 
duce a finite gluon-gluon correlation length, i.e. mag- 
netic mass. The gluon propagator has been extensively 
studied using SDE techniques |T3HIS1 and it is therefore 
worth examining to what extent topological features are 
manifested in the Green's functions. Examples of such 
studies in the context of the gluon propagator, can be 
found in ^7M\. 

In this work perform such a study in simple models 
where SDE solutions can be compared to exact results. 



In particular with these models we will be able to address 
the adequacy of truncated SDE in capturing the under- 
lying, nontrivial properties of the vacuum. In more real- 
istic models with nontrivial topology i.e. the Schwinger 
model [3TJ[55] or the abelian Higgs model Green's 
functions have been studied using semiclassical approx- 
imations by introducing dual variables that account for 
the topological defects. Here instead we will use models 
in which we can compare SD, semiclassical and exact re- 
sults. In particular we consider the following three mod- 
els that we design to capture, in a much simplified way, 
some characteristic properties of a topological vacuum of 
the more sophisticated modes, like the ones mentioned 
above. In Sec.|TT]wc compare truncated SDE results with 
the exact solution of D = + dimensions models that 
have either a unique vacuum or degenerate vacua. In 
Sec. |III| we consider a model with a quasi-periodic vac- 



uum and finally in Sec. IV we discuss the role of boundary 



conditions following the example of a particle on a circle 
i.e. a D = + 1 field theory. A D = + 1 dimensional 
theory describes a quantum particle at finite temperature 
or equivalently classical statistical mechanics of a string. 
AD = 0+0 "theory" may be considered as dimensionally 
reduced, heavy mass limit of a D = -I- 1 model. In the 
following, however, we will focus on comparing results of 
various approximation schemes, including SDE's rather 
then on their physical interpretation. Conclusions and 
outlook are summarized in Sec. |Vl 



II. UNIQUE VS DEGENERATE VACUUM 

In _D — -f dimensions the generating functional 
becomes a function of a single source variable, j and given 
by a one-dimensional integral 



Z{j) = e^(j") = / dxe-5(^)+^^ 
For the action S{x) we take 
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and depending on e consider both unique and degener- 
ate vacua: if e — +1 the action has a single minimum 
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while for e = — 1 there are two degenerate minima at 
X = iby^G/A with the a; = point corresponding to a 
local maximum. In 13 = + 1 euclidean dimensions, 
X ~> x{t) and the action dT[x'^{T)/2m+S{x{T))] (with 
S{x) defined by ([2])) describes thermal fluctuations of a 
quantum particle, which in the m — > cx) limit reduces 
to considerations of integrals as the one given in Eq. ([ij . 
We have assumed that the integral over x is not restricted 
i.e. X runs over the interval (— cx),+cx)). The semiclas- 
sical approximation will be valid in the limit of small 
coupling A — > 0, as can be easily seen once x is rescaled 
via, a: — > x = X^^^x, 



^, , I f x^ x^ 



(3) 



In the small-A limit the integral should be well approxi- 
mated by the contributions from the saddle points. For 
e = +1 there is one saddle point at x = 0, which in 
higher dimensions corresponds to a unique vacuum. For 
this reason, in the following we refer to saddle points as 
vacuum contributions, which can be unique, for e = +1 
or multiple as in the case of e = — 1 (and in the more gen- 
eral case considered in Sec. Ill I. For e = +1 the saddle 



point approximation is equivalent to the leading order 
standard, perturbative expansion in powers of A. In par- 
ticular, the two point correlation 



ct' log Zij) 



df 



(4) 



can be easily computed by expanding Eq. ([iJ in powers 
of A with the result 



.^)^i4a+?a^-^a3+?a^-^^^^ 
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A^ + 0(A''). (5) 



Comparison between the exact, numerical evaluation of 
(x^) and the above perturbative series is shown in Fig. l] 
As expected, as A decreases, the accuracy of the saddle 
point approximation improves, however, even with cor- 
rections up to O(A^) the perturbative expansion is ac- 
curate only for very small couplings, A < 0.25. This is 
an indication of the non- analytical behavior of Z{j) at 
A = 0. A similar behavior is also expected in QCD. For 
larger values of the coupling any reasonable approxima- 
tion must therefore, at least partially, re-sum the pertur- 
bative series to all orders. Since at any order of trunca- 
tion in the number of effective interactions, Schwinger- 
Dyson equations do sum up an infinite number of inser- 
tions of A one expects that a solution of a truncated set 
of SDE's will be a better approximation compared to the 
truncated perturbative expansion of Eq. ([s]). The SD 
equations follow from the identity 



F' 



1 



y 



= -r^'Hy) 



r(2)(y) dy 

where T{y) is the effective action defined by 

dW{j) 



r(y) = W{j) - jy, with y 



dj 



(6) 
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With W{j) given by Eq. Eq. ^ leads to the master 
equation, 

A / r(3)(zj) 3y 



£2/ 



6 V[r(2)(y)]3 r(2)(y) 



-y')=-rW(2/) 



(8) 

from which expectation values of any function of x, can 
be generated by taking appropriate number of deriva- 
tives. In particular the SDE for the two-point correla- 
tion. 



,2, ^ d^W{j),=o ^ 1 

^ ^ df r^'^ 
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is obtained by taking the first derivative of Eq. ([8| and 
setting the source term to zero. This gives 



(10) 



Since j — implies y = all terms odd in y, in the effec- 
tive action, Eq. ([t]) vanish and the SD equation for F'^'*) 
is obtained from Eq. Q by taking two more derivatives 
of the master equation. 



-.(4) 
- 

A 



Here Fq"-* = d"'T{y)y=Q/dy^^ is the dimensionless coupling 
in the n-point vertex of the effective action. Similarly, 
one can derive equations for all higher order vertices, 
e.g Fq"-*, n > 6 by taking more derivatives of the mas- 
ter equation. Most truncation schemes in applications of 
SDE's are based on neglecting all but a lowest few ver- 
tices. The lowest order (LO) approximation is obtained 
by setting Fg'*'' = —A i.e. neglecting dressing of the bare 
vertex implied by the rhs of Eq. (11) as well all higher 



(11) 



order vertices since they are generated by higher order 
loops. In the next to leading order (NLO) one would 
keep loop dres sing of the bare vertex, appearing on the 
rhs of Eq. ( |ll| as terms containing Fg*^ , while continu- 
ing to neglect vertices generated by higher order loops, 
Fg"-* = 0, n > 6. In LO the two-point correlator is there- 
fore given by a solution of, 



e+lix^) 



fix')'' 



(12) 



while in the NLO one needs to solve a set of the two 
coupled non-linear, algebraic equations, Eq. ( 10 ) and 



Eq. ([ll| with f[)' 



(6) 



0. It can be easily verified that 
the LO solution has the following expansion in powers of 
the coupling constant A 



c2) = l-iA + ?A2--A3 



(13) 
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i.e. it agrees with the exact result only up to the second 
order, while for the NLO solution one finds, 

599, 
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A=' + 0(A'') (14) 
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which agrees with Eq. ([5]) up to the fourth order. In 
Fig. [T] we also compare the saddle-point approximation 
of Eq. ([5| with the LO and NLO solution of SDE's for 
e = +1. Clearly the SDE equations result in the two- 
point correlation that is significantly more accurate then 
the saddle point approximation (aka perturbation the- 
ory) for large values of A. It is worth noting that this not 
because the solution of the SDE's and the perturbative 
series match to some high order in A; clearly the per- 
turbative expansion is quite inaccurate except for very 
small A. The agreement between solutions of SDE's and 
the exact result originates from the effective all order 
re-summation of the perturbative series generated by the 
nonlinear SDE's. Indeed comparing the coefficients of the 
first few terms in perturbative expansion beyond the or- 
der where SED's match the exact expansion, e.g. 619/48 
vs 599/48 one finds the difference to be only a few per- 
cent. 

For strong coupling, A > 1, however, the SDE's trun- 
cated so that higher vertices (loops) are removed, fail. 
This is because higher order vertices in the effective ac- 
tion, Fq"^ grow with A. The perturbative expansion near 
the X — Q saddle point can nevertheless be set up by 
expanding the integrand in Eq. ([T]) in terms of the "ki- 
netic term", ex^/2, instead of the interaction term, i.e. 
by expanding in powers of 1/\/A. This leads to 
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FIG. 1. Comparison of perturbative (dashed lines), 
Schwinger-Dyson (dashed-dotted line) and exact, numerical 
evaluation of the two-point correlation (solid line) for the ac- 
tion with a unique classical vacuum (e = 1). On this scale the 
NLO SDE solution is indistinguishable from the exact result. 



where all higher order coefficients, c„ 
analytically. It then follows from Eqs. 
large- A 



can be computed 
(Tol), dnl) that for 



-(4) 

- 

A 



-0.2960 



p(6) 
^ 

A3/2 



1 



-0.6276 + Ci( 

VA' 

(16) 

i.e. there is no suppression of higher order vertices, 
F(")/A"/'' = 0(1). For e = -1, SD equations truncated 
at any finite order cannot reproduce the exact result. 
This is because vertices in the SDE's generated from the 
master equation originate from expansion around x = 0, 
which is a local maximum and a metastable sate in higher 
dimensions. For e = — 1, {x^) is still positive and as a 
function of A it is non- analytical at A = 0, where it has 
a pole, while the LO solution Eq. ( 12 ) is analytical and 
for A = gives {x'^)lo = — l! The action has two min- 
ima and in the saddle point approximation the integral in 
Eq. ([T]) is approximated by a sum of gaussian fluctuations 
around each of them with the difference between the full 
action and gaussian approximation treated as perturba- 
tion. This leads to 

(17) 

and is compared to the exact, numerical result in Fig. [2] 




FIG. 2. Saddle-point approximation computed to a varying 
order in A (dashed lines) compared with the exact result for 
e = 1 obtained by numerical integration (solid line). The 
□ (A*) perturbative results is indistinguishable form the exact 
result. For easier comparison we subtracted the leading, 6/A 
term from {x^). On this scale the solution of SD equation 
tends to —6/A — —oo as A — > 0. 
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III. MULTIPLE QUASI-DEGENERATE VACUA 



In QCD large gauge transformations are not con- 
strained by the Gauss' law and result in topologically 
disconnected field configurations I27ti31j . In the semi- 
classical approximation these configurations correspond 
to degenerate classical vacua and the true vacuum state 
is a linear superposition of these vacua. To access the 
applicability of SDE equations in the case of such multi- 
ple saddle points, we again make a simple model for the 
action in Z? = -f 0. Since there is no tunneling without 
the "time" direction, it is possible that SDE's equations 
in the realistic case with instantons are more accurate 
then the one ones in the no-time model discussed here. 

To mimic the effect of multiple vacua in the partition 
function we consider the following action 
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S{x) = — 
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2(1 - cos(a;)) 



X 

A2 



(18) 



For simplicity, to be able to compare with SDE's we de- 
fine X in the (— oo,-|-oo) interval. The multiple saddle 
points originate from the cos term the role of the second 
term on rhs in Eq. 18 is to make the integral in Eq. ([T]) 
well defined. Because of this term the minima of S(x) 
are not really degenerate, and their relative contribution 
depends on A. For weak coupling, g < 1 the expectation 
value of (x^) is determined in the saddle point approxi- 
mation by the number of minima of the action. For fixed, 
finite A, (x^) — > when g — since in this case the term 
x"^ /{g^A"^) is large for all minima of cosx except the one 
at a; = 0. In this limit, therefore, the integral becomes 
dominated by the single minimum at a; = and one finds. 



2„2 



A'g 



2(A2+52) 



2 



(19) 



Since the problem becomes effectively that of a single 
vacuum one expects the SDE to yield a similar result. 
This indeed is the case as shown in Fig. |4j As A increases 
for fixed 5 < 1 the contribution from the minima of cos x 
at X — 2'Kn, \n\ > are no longer suppressed and in the 
limit A — > 00 one obtains 



A2 
2 ■ 



(20) 



This is an interesting limit, since even though there are 
multiple vacua contributing to the partition function in- 
tegral their contribution is approximately equal to that 
of a broad minimum given by x^ j {c.f. Fig. p3|, a nd it 
is this gaussian distribution that results in Eq. (|20p. 

The SD equations for the action of Eq. ( 18 1 are derived 
from the operator identity 



sm 



9^ 



1 d 

TM(y)dy 



+ % = -^''^iy)- (21) 



In particular for the two-point correlation in the LO ap- 



FIG. 3. A typical action used here in a model study of 
multiple, quasi-degenerate vacua. For large A and/or small 
g the the number of minima contributing to the partition 
function is proportional to h? / . For the case shown in this 
figure, besides the central minimum there are additional A'^ = 
32 local minima. 



proximation, T^'^'^{y) = Fq^-*, it yields 



2 d 
^2 dy ' 



r(j2) dy 



A - _r(2) 
A2 - ^0 



(22) 



The non-linear, algebraic equation for {x'^) — — 1/Fg 
can be derived from the above by noticing that the argu- 
ment of the sine can be written as an expectation value 
of ladder operators 



dy 



sm 



1 



p(2) dy 



v/(^(0|asin 



with a = 1/ yjjx^d/dy and 
to the following gap equation 



2 _i£!i 



2 

A2 



|0) 
(23) 

^ (x2)y. It finally leads 



(24) 



It can be verified that for finite A in the g — >■ limit 
the solution of Eq. ( 24 1 agrees with Eq. ( 19 1 while in the 
limit A — cx) it agrees with Eq. ( 20 1 In the former case, as 
discussed above, the SD equation reproduces the result 
of the problem with a unique vacuum. For large-A on 
the other hand, it indeed follows from Eq. (24) that even 



though multiple saddle points contribute the net effect 
is equivalent to that of a single minimum with a large 
width ^ A leading to solution for {x^) that agrees with 
Eq. ( 20 ) , as shown in Fig. |4] 

Outside of the two limits, several local minima of he 
action contribute to the partition function and the LO 
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FIG. 4. Comparison 
Schwinger-Dyson 



Eq. 



between exact (solid line) and 



( 24 1 appr oxi mation (dashed line) to 
(a;^) for the action given by Eq. ([l8|. We used = 0.01 

(c/. Fig|3|. The dashed-dotted 
limits given by Eqs. ([T9|, ([20| 



ines correspond to the two 
were, as discussed in the 
text Schwinger-Dyson approximation becomes exact. For the 
rang e o f couplings shown the saddle point approximation of 
Eq. ( 25 I is indistinguishable from the exact result and involves 



summation over A'' = 33 saddle points. 



SD equation fails. This is seen in Fig. |4] for g in the mid- 
range, 0.2 < 5 < 1. In this range, where multiple minima 
contribute one expects summing over integrals around all 
saddle points is a better approximation to the generating 
functional compared to the SDE. In this, saddle point 
approximation 



(25) 



and comparison with the exact, numerical result is shown 
in Fig. [sj Indeed the saddle point approximation works 
while SDE fails in this range of couplings. 



IV. PARTICLE ON A CIRCLE 

In QCD the domain of gauge fixed gluon field, i.e. 
the fundamental modular region may be non-flat with 
a non-trivial measure specified by the Fadeev-Popov de- 
terminant [32 . This highly complicates the Schwinger- 
Dyson formalism, which needs to take into account the 
boundary of fundamental modular region or the Gribov 
horizon |33| . In the following example we investigate 
what happens in a model in which the dynamical vari- 
able has a non-trivial boundary under an approximation 
when this boundary is ignored. For this purpose we con- 
sider quantum mechanics of a particle on a circle |34j i. e. 
D = Q + 1 dimensional field theory. The variable x{t) 



FIG. 5. Comparison between exact, numerical computation 
and the saddle point approximation to {x^) for different val- 
ues of / h.^ . As A increases for fixed g so does the num- 
ber of saddle points. The three cases shown correspond to 
A'' = 3,7,33 points respectively (c/. Eq. ([25|). As A in- 
creases, saddle points become degenerate and Eq. |25| becomes 
an increasingly better approximation. 



now describes location of the particle on a unit circle 
Q < x{t) < 27r as a function of time and for the Hamil- 
tonian we choose the free kinetic energy for a particle of 
unit mass 2H — = —d^. Since the manifold is com- 
pact, the wave function must satisfy the boundary condi- 
tion ip{x) = exjp{i9)ip(x+2TT) and in the following we take 
9^0. The normalized eigenvectors of the Hamiltonian 
are spanned by ipmix) = exp{imx) / with integer m 
and the corresponding energies are Em = rr? 12. The vac- 
uum expectation value at the euclidean time (t = —il3), 
i.e. temperature-dependent correlation function is then 
given by 



{x{(3)xi0)) = (0|xe-^^a;|0) = + 



(26) 



At low temperatures, T = 1//3 — the correlation func- 
tion is dominated by the lowest energy quantum sate and 
{x{oo)x{0)) — ?> TT^. In this case the restriction that x be 
on a circle is important. At high temperatures, however, 
the system becomes semiclassical, and the particulars of 
the topology of the quantum system should become ir- 
relevant. Also in this limit expectation values should be 
well approximated by contributions from small amplitude 
fluctuations around solutions of the classical equation of 
motion. In this case truncated SD equations should also 
be a good approximation. In our simple example we as- 
sume no interaction thus the SD equation and the semi- 
classical, saddle point approximation give the same re- 
sults and both of them pertain to a formulation of the 
problem in terms of the variable dual to quantum number 
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m [3S]. This variable is just is the classical coordinate x 
and the duality transformation m o a; is given by 



dx 



' 2/3 " 



(27) 



leading to 



dx' dx 

{x{l3)x{0)}= ^ :z:'(x',/3|x,0)x (28) 
Jq Ztt Ztt 



where 



{x ,/3|a;,0) = a/ — 2^ exp ' 



2/3 



(29) 



At finite temperature, q counts the number of times the 
particle wraps around the circle. The duality between m 
and x is clear; at high temperature, /? — ?> and x is well 
defined, (while m is not) since 



(a;',/3|x,0) ^ 2t:5{x' - x) 



(30) 



we have done such comparison in simple models to illus- 
trate strengths and potential limitations of the SDE ap- 
proach. In particular we have shown that if the partition 
function integral is dominated by a single semiclassical 
configuration, (one saddle point of the action), even the 
one-loop SD equation gives a much better approxima- 
tion for a two-point correlation than perturbative series 
truncated at some high order. In the case of multiple- 
saddle points, however, we have shown that the semi- 
classical approximation, when applicable is more reliable 
while the SDE approach may not necessarily be captur- 
ing the correct physics. In such cases SDE's can, however 
be successful when formulated in dual variables. Clearly 
the models studied here are very naive and more realistic 
problems need to be considered before definite conclu- 
sions about applicability of truncated SD equations to 
describe physics of confinement are made. 



with only the term with q ~ Q contributing. It then im- 
mediately follows that, (a;(0)a;(0)) = 47r^/3. In the high 
temperature limit, the quantum variable m is not well 
defined, i.e. in the sum in Eq. (261 an infinite number 
of terms contribute to give 47r^/3. At low temperatures, 
on the other hand, the system is quantum. As /3 — > cx), 
m becomes well defined, to (while x is not) and 
(a;(oo)a;(0)) — ^ tt^. To obtain the low temperature limit 
of the correlation function using the classical represen- 
tation of Eq. (291 it is necessary to integrate over the 
entire range of x and x' . Furthermore it is necessary to 
allow X (or x') to wrap around the circle an arbitrary 
number of times i.e. sum over all q ^ N. The approxi- 
mation in which the topology of the boundary is ignored 
corresponds to retaining only the q — term in the sum 
in Eq. (29), which, as follows from the discussion above. 



should work fine in the high temperature limit but fail at 
low temperatures. This is shown in Fig. |6] 
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FIG. 6. Comparison between exact (dashed line) and ap- 
proximate evaluation of the correlation function as a function 
of temperature. In the approximation Eq. 1 29 1 only the q = 
term is retained. 



Phenomenological applications of Schwinger-Dyson 
equations require truncations in the number of retained 
effective interactions. It is important to access the reli- 
ability of such truncation in QCD where Green's func- 
tions are not necessarily dominated by small fluctuations 
around the perturbative vacuum. In particular, confine- 
ment is expected to be related to topologically non-trivial 
field configuraiton like magnetic monoples and vortices 
that lie on the Gribov horizon [37^ and produce tunnel- 
ing between the degenerate, local minima of the action. If 
the effective potential for gauge degrees of freedom were 
known in principle one could compare the SDE approach 
with other, i.e. semiclassical approximations, (as can be 
done for example in the case of compact QED PB]). Here 



7 



ACKNOWLEDGMENTS FG0287ER40365 and by DFG under contract DFG- 

Re856/6-3. 

This work was supported in part by the US 
Department of Energy grant under contract DE- 



[1] C. D. Roberts and A. G. Williams, Prog. Part. Nucl. 

Phys. 33, 477 (1994). 
[2] R. Alkofer and L. von Smekal, Phys. Rept. 353, 281 

(2001). 

[3] P. Maris and C. D. Roberts, Int. J. Mod. Phys. E 12, 
297 (2003). 

[4] R. Alkofer, C. S. Fischer and F. J. Llanes-Estrada, Phys. 

Lett. B 611, 279 (2005) [Erratum-ibid. 670, 460 (2009)]. 
[5] R. Alkofer, C. S. Fischer, F. J. Llanes-Estrada and 

K. Schwenzer, Annals Phys. 324, 106 (2009). 
[6] L. Del Debbio, M. Faber, J. Greensite and S. Olejnik, 

Phys. Rev. D 55, 2298 (1997). 
[7] K. Langfeld, H. Reinhardt and O. Tennert, Phys. Lett. 

B 419, 317 (1998). 
[8] M. Engelhardt, K. Langfeld, H. Reinhardt and O. Ten- 
nert, Phys. Rev. D 61, 054504 (2000). 
[9] J. Greensite, Prog. Part. Nucl. Phys. 51, 1 (2003). 
[10] Y. Nambu, Phys. Rev. D 10, 4262 (1974). 
[11] S. Mandelstam, Phys. Lett. B 53, 476 (1975). 
[12] A. M. Polyakov, Nucl. Phys. B 120, 429 (1977). 
[13] G. 't Hooft, Nucl. Phys. B 190, 455 (1981). 
[14] L. von Smekal, A. Hauck and R. Alkofer, Annals Phys. 

267, 1 (1998) [Erratum-ibid. 269, 182 (1998)]. 
[15] C. S. Fischer, A. Maas and J. M. Pawlowski, Annals 

Phys. 324, 2408 (2009). 
[16] A. C. Aguilar, D. Binosi and J. Papavassiliou, Phys. Rev. 

D 78, 025010 (2008). 
[17] J. M. Cornwall, Phys. Rev. D 57, 7589 (1998). 
[18] F. V. Gubarev, M. I. Polikarpov and V. I. Zakharov, 

Phys. Lett. B 438, 147 (1998). 
[19] M. N. Chernodub, M. I. Polikarpov and V. I. Zakharov, 
Phys. Lett. B 457, 147 (1999). 



[20] A. P. Szczepaniak and H. H. Matevosyan, Phys. Rev. D 

81, 094007 (2010). 
[21] C. Adam, Czech. J. Phys. 46, 893 (1996). 
[22] T. Radozycki, Phys. Rev. D 60, 105027 (1999). 
[23] P. Orland, Nucl. Phys. B 428, 221 (1994). 
[24] M. Sato and S. Yahikozawa, Nucl. Phys. B 436, 100 

(1995). 

[25] E. T. Akhmedov, M. N. Chernodub, M. I. Polikarpov and 
M. A. Zubkov, Phys. Rev. D 53, 2087 (1996). 

[26] M. I. Polikarpov, U. J. Wiese and M. A. Zubkov, Phys. 
Lett. B 309, 133 (1993). 

[27] R. Jackiw and C. Rebbi, Phys. Rev. Lett. 37, 172 
(1976). 

[28] A. Belavin, A. Polyakov, A. Schwartz, and Yu. Tyupkin, 

Phys. Lett. 59, 85 (1975). 
[29] H. Reinhardt, Nucl. Phys. B 503 (1997) 505. 
[30] C. Ford, U. G. Mitreuter, T. Tok, A. Wipf, 

J. M. Pawlowski, Annals Phys. 269, 26-50 (1998). 
[31] O. Jahn, F. Lenz, Phys. Rev. D58, 085006 (1998). 
[32] V. N. Gribov, Nucl. Phys. B 139, 1 (1978). 
[33] D. Zwanziger, Nucl. Phys. B 323, 513 (1989). 
[34] N. Fjeldso, J, Midtdal, and F. Ravndal, J. Phys. A: 

Math. Gen. 21, 1633 (1988). 
[35] T. Banks, R. Myerson and J. B. Kogut, Nucl. Phys. B 

129, 493 (1977). 
[36] S. D. Drell, H. R. Quirm, B. Svetitsky and M. Weinstein, 

Phys. Rev. D 19, 619 (1979). 
[37] J. Greensite, S. Olejnik, D. Zwanziger, Phys. Rev. D69, 

074506 (2004). 



